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One way of deriving the atmospheric extinction coefficient from lidar measure- 
ments is to start from the single scattering lidar equation, which we will write 
in the form: 

(1) P(r)=P()A — Mill exp (’-/!' ^(r^dr 1 ] 

2 r 2 0 

This equation cannot be solved directly, since both e(r) and <r(r) are unknown. 
In the approach under consideration here, one attempts to solve the equation by 
postulating a relation between g and a of the form: 

(2) B(r)=c [<r (r) ] k 


Such a relation cannot hold in general, but may be approximately valid for a 
wide range of realistic atmospheric conditions. Real atmospheres will deviate 
from this relation, of course, so it is important to know how large are the 
errors which we should expect due to such deviations. 

In the present study this question has been addressed by conducting simulations 
of lidar experiments and extinction coefficient calculations and examining the 
resulting errors. In order to make the simulations realistic, the volume back- 
scatter coefficients and extinction coefficients used in the simulations were 
based on measured particle distributions, in this case on the particle size dis- 
tributions measured in Meppen, Germany in the fall of 1980 by j. D. Lindberg 
and others of the Atmospheric Sciences Laboratory (Lindberg, 1982). 

In that experiment particle size distributions were measured with a balloon- 
borne particle spectrometer (a PMS FSSP 100) during balloon descents from 
about 600 meters above ground to the surface. Each particle spectrum represents 
a layer of at most a few meters thickness so an altitude resolved profile of 
particle size spectra is measured. 

T. L. Barber (Barber, 1985) used a Mie code to compute linear extinction co- 
efficient and volume backscatter coefficient profiles from the particle size 
spectra for each balloon sounding. In these computations it was assumed that 
the individual droplets consisted purely of water. 

Simulated lidar returns were constructed by inserting the computed values of 
extinction and backscatter coefficients into a discrete version of equation 
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(1). In order to mimic the range gating characteristics of an actual lidar 
used by the author (Measure and Rubio, 1982) at the Meppen experiment, the 
interval in the discrete version of the lidar equation was chosen to be 
three meters, and the backscatter and extinction coefficient profiles were 
interpolated at the chosen points. This resulted in a set of values P(r) 
for r = 3n meters, n = 1, 2, . . . which constituted the simulated 
lidar return. (For simplicity, the factor PqCXA/ 2 was set to 1). 

This simulated lidar return was then used to compute the extinction coeffi- 
cient by Klett's algorithm (Klett, 1981): 

(3) c(r ). e*p[(S(r)-S(r m ))A] 

[°m (''m*]" 1 ex PC<S(r)-S(r m ) )/k] 

k Jr 

2 

where S(r) = ln(r P(r)). This algorithm requires that the boundary value 
o(r ) be specified. In the simulation, the value used was that computed 

from the particle spectra, so that the simulations were "perfect" at the 
point r^ and the only error should be that due to the deviation from the 

backscatter to extinction relation (2). 

Figure 1 indicates the data flow in the simulation. Representative results 
of the simulation are indicated in figures 2-4. Note that using a k value 
other than 1 may result in a significantly smaller error. Note also., that 
no other potential sources of error are considered in the simulation. 
Absorption, multiple scattering, and instrumental errors all may contribute 
to error in a real lidar measurement but are not represented here. 
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